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Abstract 

We prove the regularity for weak solutions to a class of ul- 
traparabolic equation, with measurable coefficients. The results gen- 
eralized our recent regularity results of Prandtl's system to high 
dimensional cases. 

keywords: Ultraparabolic equations, Moser iteration, regularity 

1 Introduction 



The ultraparabolic equation arises in many applications, for example, 
fluid dynamics, mathematical finance, degenerated diffusion process, etc. 
There are more and more studies on this problem in recent years. The reg- 
ularity of this type of equation becomes interesting since it has some special 
algebraic structures and is degenerated. It is still unclear in general, whether 
the interior C° regularity results hold for weak solutions of the ultraparabolic 
equations with bounded measurable coefficients like the parabolic cases. 



*The author currently is working at NSFC 
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In the study of boundary layer problem, we obtained the existence global 
weak solution in the class that Oleinik considered under the assumption that 
the pressure is favorable [13]. One of the interesting question is whether 
the weak solution is actually smooth? In particular, in the two dimensional 
Prandtl's system with constant pressure, in the Crocco variable, we can de- 
duce the following equation 

du du 2'^^ g 
dt ^ dx dy^ 

This is of strong degenerated parabolic type equations, more precisely, an ul- 
traparabolic type equation. However, it satisfies the well known Hormander's 
hypoelliticity conditions, which sheds lights on the smoothness of weak so- 
lutions. It is proved by Polidoro and Ragusa [12] that the weak solution is 
in the class, if the coefficient is in the class of VMO. That is, if the weak 
solution u is also in the class of VMO, then the solution u is continuous 
and then the result of Xu [15], or Bramanti, Cerutti and Manfredini [1] im- 
plies that the solution is also smooth. In our case, however, we can only 
prove that the weak solution is in the class of BV space. It is interesting 
if the weak solution of equation (1.1) is still smooth when the coefficient is 
only measurable functions. 

On the other hand, equation (1.1) has the divergent form if we replace u 
by A recent paper by Pascucci and Pohdoro [11] proved that the Moser 
iterative method still works for a lass of ultraparabolic equations with mea- 
surable coefficients including equation (1.1). Their results showed that for a 
non- negative sub-solution u of (1.1), the norm of u is bounded by the LP 
norm of w (p > 1). This is a very important step toward the final solution of 
regularity of the ultraparabolic equations. 

We proved in [14] that the weak solution that we obtained in [13] of (1.1) 
is of class and then u is smooth. In this paper, we are concerned with the 
regularity of solutions of the ultraparabolic equations. We shall generalize 
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the result in [13] to high dimensional cases in this paper. 

We consider a class of Komogorov-Fokker-Planck type operator on 

mo N 

(1.2) Lu= d:cj {aij{t, x)d:ci u)+ Y bijXid^. u-dfU^O, 

i,j=l i,j=l 

where (x^t) e 1 < mo < -/V, and bij is constant for every i,j — 

1, - • • ,N. We make the following assumptions on the coefficients of L: 

{Hi) = ttji e L°°{R^~^^) and there exists a A > such that 

1 mo mo mo 

i=l i,j=l i=l 

for every {t, x) G and { G 

{H2) The matrix B — (6ij)ivxiv has the form 
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where B^. is a matrix ruk-i x with rank rrih and mo > mi > • • • > m^, 
mo + mi + • • • + nid = N. | | < A where the norm 1 1 ■ 1 1 is in the sense of 
matrix norm. 

The requirements of matrix B in {H2) ensures that the operator L with 
constant aij satisfies the well-known Homander's hypoellipticity condition. 

The Schauder type estimate of (1.2) has been obtained. Besides, the reg- 
ularity of weak solutions have been studied by Bramanti, Cerutti and Man- 
fredini [1], Manfredini and Polidoro [7], Polidoro and Ragusa [12] assuming 
a weak continuity on the coefficient aij. It is quite interesting whether the 
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weak solution has Holder regularity under the assumption (Hi) on aij. The 
first advances is the work of A. Pascucci and S. Polidoro [11] who proved that 
Moser iteration method still works for equation (1.2). One of the approach 
to the Holder estimates is to obtain the Harnack type inequality. In the case 
of elliptic equation with measurable coefficients, the Harnack inequality is 
obtained by J. Moser [8] via an estimate of BMO functions due to F. John 
and L. Nirenberg together with the Moser iteration method. J. Moser [9] also 
obtained the Harnack inequality for parabolic equations with measurable co- 
efficients by generalizing the John-Nirenberg estimates to the parabolic case. 
Another approach to the Holder estimates is given by S. N. Kruzhkov [6], 
[7] based on the Moser iteration to obtain a local priori estimates, which 
provides a short proof for the parabolic equations. 

We prove a Poincare type inequality for non-negative weak sub-solutions 
of (1.2). Then we apply it to obtain a local priori estimates which implies 
the Holder estimates for ultraparabolic equation (1.2). 

Let Dmo be the gradient with respect to the variables xi,X2, - ■ ■ , Xmo- And 

TV 

We say that m is a weak solution of (1.2) if it satisfies (1.2) in the distribution 
sense and u, D^^u, Yu G Lf^^. 

Our main result is the following theorem. 

Theorem 1.1 Under the assumptions (Hi) and {H2), the weak solution of 
(1.2) is Holder continuous. 
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2 Some Preliminary Results 



One of the important feature of equation (1.2) is that the fundamental solu- 
tion can be written down explicitly if the coefficients a,;,- is constant, (see [1], 
[4]). Besides, there are some geometric and algebraic structures in the space 
jl^+i induced by the constant matrix B (see for instance, [1]). 

Let E{t) — exp{—TB'^), where E{t) is a polynomial of degree d in r with 
N X N matrices coefficients. For {t,x), {T,y) e R^'^^, set 

{t,x)o{T,y) = {t + T,y + E{T)x). 

Then (i?^+^,o) is a group with neutral element (0,0); the inverse of an 
element {t,x) is {t,x)~^ = {—t,—E{—t)x). The left translation by (r, y) 
given by 

{t,x) ^ {T,y) o {t,x), 
is a invariant translation to operator L when coefficient aij is constant. 
The associated dilation to operator L with constant coefficient aij is given 

by 

6x = diag{X'^, Xlmo, A^/mi, ■ ■ • , X'^'^'^^ Ima) ^ 

where I,n^. denotes the x nik identity matrix. Then the operator is homo- 
geneous of degree 2 with respect to the dilation Sx. Let 

Q — mo + 3mi -|- • • • -|- {2d + l)md- 

Then the number Q + 2is usually called the homogeneous dimension of R^'^^ 
with respect to the dilation S\. 

The norm in related to the group of translations and dilation to 

the equation is defined by 

||(t,a;)||=r 
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if r is the unique positive solution to the equation 

— ^ H ^ H \ — H = 1 

where (i, x) e R^+^ \ {0} and 

'3^1 ■ ■ ■ Olmo 1) C^TOo+l ' ' ' OijYiQ+mi 3, • • • , 
<^mo+--+md-i+l — ■ ■ ■ — OlN — 2d + 1. 

And 11(0, 0)11 = 0. The balls at a point (to, Xq) is defined by 
Brito,Xo)^{{t,x)\ \\{t,x)-^o{to,xo)\\<r}. 

Let 

l3~{to,xo) = Br{to,xo) n{t < to}. 

For convenience, we sometimes use the cube replace the balls. The cube 
at point (0, 0) is given by 

Cr{0,0) = {{t,x)\ \t\<r^, |xi| < r"\---,|a;jv| < r"^}. 

It is easy to see that there exists a constant A such that 

Cr(0,0) CB,(0,0) cCa,(0,0), 

where A only depends on B and N. 

When the matrix (ajj)7Vx7v is of constant matrix, we denoted it by Aq. 
Then the operator L takes the form 

Lo = div{AoD) + Y. 

We let z — {t,x). The fundamental solution ro(-,C) of Lq with pole in 
C e R^'^^ has been constructed (see [1]) as follows: 

ro(^,C) = ro(C-'o^,o), ^,Cei?^+\ z^^c, 
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And ro(^,0) can be written down explicitly. There are some basic estimates 
for To 

ro(^,C)<c||r'oz||-^, 

\d,To(z,C)\<c\\C'oz\\-'^-\ 

where i — 1, - ■ ■ , mo. 

A weak sub-solution of (1.2) in a domain i7 is a function u such that u, 
DmoU, Yu e L}oSS^) and for any e (7^(0), > 0, 



(2.1) / (\)Yu-{DufAD4)>^. 

Jn 



A result of Pascucci and Polidoro obtained by using the Moser's iterative 
method (see [11]) states as follows. 

Lemma 2.1 Let u be a non-negative weak sub-solution of (1.2) in Vt. Let 
{tQ^Xo) G and B~{tQ,Xo) C and let p > 1. Then there exists a positive 
constant c which depends only on A and the homogeneous dimension Q such 
that, forO<r<l 



(2-2) sup u^<^,L ^ 

Br{t0,X0) ^ JBr(to,Xo} 
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provided that the last integral converges. 

We copy a classical potential estimates (see [3]) here to prove the Poincare 
type inequality. 

Lemma 2.2 Let a E (0,(5 + 2) andG G C(i?^"''^ \ {0}) be a 5\-homogeneous 
function of degree a — Q — 1. If f E L^{R'^'^^) for some p G (0, oo), then 



Gf{z) ^ J^^^^ GiC' ■ z)f{OdC, 



7 



is defined almost everywhere and there exists a constant C = C{Q,p) such 
that 

(2-3) \\Gf\\L,(RN+i) <C m&x \G{z)\ 1 1/| |LP(Riv+i), 

where q is defined by 

1 _ 1 _ a 
q~ p Q + 2' 

3 Proof of Main Theorem 

We want to obtain a local estimates of solutions of the equation (1.2), for 
instant, at point {to,xo)- Since the equation (1.2) is invariant under the left 
group translation when a^- is constant, we may consider the estimates at a 
ball centered at (0,0). We mainly prove the following Lemma 3.4 which is 
essential in the oscillation estimates in Kruzhkov's approaches in parabolic 
case. Then the C° regularity result follows easily by the standard arguments. 

For convenience, in the following discussion, we let ) and 

X = {x',x). We consider the estimates in the following cube, instead of , 

Cr = {{t,x)\ -r^ <t<0, \x'\ < r, \xmo+i\ < AA^V • • • , |a;jv| < AA^V^'^+^}. 

Let 

Sr^{x I Ix^o+il < AiVV3---,|a;jv| < AAtV^+^I. 

Let < a, P < 1 he constant and 

Kr = {x'\ \x'\ < r}, 

Sar = {x I |a;^o+il ^ XN'^iarf • • • , |a;jv| < XN'^ {arf'^^^} , 
Kar = {x'\ \x'\ < ar}. 
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Now for fixed t, let 

M = {{x',x) e Kpr X SfSr, U > h}. 

In the following discussions, we sometimes abuse the notations of B~ and 
Cr, since there are equivalent. 

Lemma 3.1 Suppose that u{t,x) > be a solution of equation (1.2) in B~ 
centered at (0,0) and 

mes{{t,x) & , u>l}> ^mes{B^). 

Then there exist constant a, (5 and h, < a, P,h < 1 which only depend on 
X and N such that for all t e {—ar^, 0), 

mes{J\ft} > -^meslK/jr x S^r}- 

Proof: Let 

V = ln+( ^ ^ 



u + n» 



where h is a constant < /i < 1 to be determined later. Then v at points 
where v is positive, satisfies 



mo 



(3.1) J2 {aij{t, x)d^, v) - {DvfADv + x^BDv -dtV^O. 

Let r]{x') be a smooth cut-off function so that 

r]{x') — 1, for \x'\ < f3r; 

r]{x') = 0, for \x'\ > r. 
Moreover, < 7/ < 1 and \D^,r]\ < 
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Multiplying rf^x') to (3.1) and integrating by parts 

!kp, !spr ^(^' x',x)dxdx' + ^ /^^^ Js^^ \Dmov\'^dxdx'dt 

(3.2) < ■^Qj^zp)2mes{S0r)mes{Kpr) + Jk, Jg^^ rfx'^BDvdxdx'dt 

+ ixr Isgr ^('^^ x',x)dxdx', 
where C only depends on A and N. 

Integrating by parts, we have for any i,j 

(3.3) / / rfxihijd^vdxdx'<^^j—,(3~'^^^^^ 

JKr JSfir 4iV lis 

Then 

(3.4) f jj^ jg rfx^BDvdxdx'dt<^(5-^^\n{^)mes{Spr)mes{Kpr). 
We shall estimate the measure of the set N't- Let 

/u(i) = mes{{x' ,x)\ x' e Kj., x & Sr, u> 1}. 

By our assumption, (for convenience, we may let B~ be replaced by C~), for 
< q; < I 



2 

That is 



-r^mes{Sr)mes{Kr) < I ii{t)dt — I it{t)dt + / it{t)dt. 

Z J—r^ J—r^ J—ar'^ 



IJi{t)dt > (- — a)r^mes{Sr)mes{Kr). 
Then there exists a r e (— r^, — ar^), such that 
(3.6) /u(t) > (- — q;)(1 — Q;)~^mes(5'r)mes(K^). 



10 



From (3.2) and (3.6), we have by noticing v = when u > 1 

(3.7) f f v{t,x' ,x)dxdx' < -{1 — a)~^'rnes{Sr)mes{Kr)ln{—). 

J Kr JS/Sr 2 }ia 

Now we choose a (near zero) and /3 (near one), so that 

By (3.2), ( 3.4), (3.7) and (3.8), we deduce 

ln{^)mes{Kfir x V \M) < [C(l - /3)-^/3-«+ 

(3.9) 

+Hn{-^)]mes{Kpr x 5^^). 

Since 

Hh-i) ^ 9 

ln(/.-i) ^8' ^« ^^0' 
then there exists constant hi such that for < /i < /ii and t e [— ar^, 0] 

mes{Kpr X SprXMt) < ^mes(X^r x S'/jr)- 
Then we proved our lemma. 

Corollary 3.1 Under the assumptions of Lemma 3.1, we can choose 6, < 
9 < a and 6 < 13 small enough so that 

mes{B-p,\Bs^ n {{t,x)\ u>h}}> Co{a,(5,K)mes{Bp,}, 

where < Co(q;,/3,A) < 1. 

Let x(s) be a smooth function given by 

x(s) = 1 if s< Ver, 
X{s) = if s> pr, 



11 



where \fd < | is a constant. Moreover, we assume that 
and < if < s < ^r. We set 

N 2 
i=mo 

Mx) = x{o\x'\), 

(3.10) (l){t,x)^Mt,x)Mx), 
where ci > 1 is chosen so that 

for -r^ <t<0 and x e Kj. x Sr- 

We now have the following Poincare's type inequality. 

Lemma 3.2 Let w be a non-negative weak sub-solution of (1-2) in Bi. Then 
there exists a constant C, only depends on A and N , such that for r < 9 < 1 

(3.11) 4(-W-^»)^-SCj^/ 

where Iq is given by 

(3.12) lo ^ max^Jhiz) + C2{z)], 
and 

(3.13) h{z)= [ [{(l),AoD(l)o,Dro{z,-))w-ro{z,-)wY(j)]{C)dC, 

C2{z)^ [ [{(t>oAoD4>^,D^o{z,■))w]{OdC, 
where Vq is the fundamental solution, and 4> is given by (3.10). 
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Proof: We represent w in terms of the fundamental solution of Tq. For 
z e Bqj., we have 

H^) = Ib- [{AoD{w(t>),Dro{z, •)) - Toiz, ■)Y{wcf>)]{C)dC 

(3.14) 

^h{z) + l2{z) + Is{z) + C2{z), 
where Ii{z) is given by (3.13) and 

h{z)= f [{{Ao-A)Dw,Dro{z,-))<p-ro{z,-){ADw,D<f>)]{C)dC, 



Is{z) = / [{ADw, D{To{z, ■)(!>)) - To{z, ■)<pYw]{C)dC: 

JBr 
F 

And C2{z) denotes the remaining parts of the integral of Ii. 



From our assumption that w isa. weak sub-solution of (1.2), then l3{z) < 
(see[ll]). Then in Bq^., 

0<{w{z)-Io)+<l2{z). 

By Lemma 2.2 we have 

(3-15) \MmB-^) < ^r||/2||^,+ 4^^_^ < -—^\\D^^w\\^,^^-y 

Then we proved our lemma. 

Now we apply Lemma 3.2. to the function 

w — Iri^ 9". 

u + hs 

We estimate the value of /q given by (3.12) and (3.13) in Lemma 3.2. 

Lemma 3.3 Under the assumptions of Lemma 3.2, there exist constants Aq, 
Tq and rQ < 6 only depend on constants a, (3, A, X, (j) and 6, < Xo < 1, 
such that for r < Tq 

(3.16) < AoMA)- 

hs 
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Proof: We first note that the support of the function ^o-D^i is contained in 
the set B~ fl {{t,x',x)\ < \x'\ < |}. Then it is easy to check that there 
exists a constant C which only depends on A and A*" such that, for z e 

\C2{z)\ < C9ln{\). 

ha 

Therefore C2{z) ^ as 9 ^ 0. 

Now we let w = 1 then (3.14) gives, for z e Bg^., 

(3.17) 1= / [{<P,AoD(l)o,Dro{z,-))-ro{z,-)Y4>]iOd(: + C2iz)\^=i, 

JB- 

where is given by (3.10). Since the matrix C~^{t) is positive definite for 
t > and by the assumption of matrix B, one can check 

Y{C-\\t\)x,x) = -{AoC-\\t\)x,C-\\t\)x), 

(3.18) y(C-^(|t|)e*^^a;,e*^^x) = -(Aoe*^C-^(|t|)e*^^a;, e*^C-^(|t|)e*^^x), 
then by the choosing of Ci, it is easy to see that 

(3.19) Y(j) < 0. 

For 2; = 0, by our construction of 0, we have 

{(l>iAoD4>o,Dro{z,-))>0, 

therefore 

(3.20) {<f>iAoD(f)o, DFoiz, •)) - ro(;^, ■)¥(!> > 0. 
In fact, 

(AoD(C-i(|t|)e*^^a;,e*^^x),D(C-^(|t|)e*^^x,e*^^a;)) > 0. 
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Wc note that the support of x'{^) is in the region v^r < s < (3r. Thus for 
some [3' < f3, the set B^,^. \ with |t| > is contained in the support 
of (/)' and then the inequahty holds in (3.20). By the choosing of Ci, we know 
that (3.20) is positive in B^,^ \ B^^ with \t\ > Then the integral of 

(3.20) on the domain B^,j.\B^^ with \t\ > is lower bounded by a positive 
constant which independent of small r and 6. 

Then we may choose 9 small so that (3.20) still holds for z e B'^^ and 
the inequality holds in B^,^ \ with \t\ > 9'^r^. Therefore the integral 
function in (3.17) is nonnegative for z e and positive in B^,^ \ B^^ with 
\t\ > 9'^r'^. Since w — when u > h and w < ln{^), then by choosing 9 
small enough, our lemma follows from Corollary 3.1 and (3.17). 

Lemma 3.4 Suppose that u{t,x) > be a solution of equation (1-2) in B~ 
centered at (0, 0) and 

mes{{t,x) G u > 1} > -mes{B~). 

Then there exist constant 9 and Hq, < 9,ho < 1 which only depend on X, 
Xq and N such that 

u{t, x) > ho in Bq^. 

Proof: We consider 

w = ln+( -j), 

u + hs 

for t e [—ar^, 0]. By applying Lemma 3.2 to w and as in the proof of Lemma 
3.1, we have 

(3.21) -f_{w-Io)l<C{ln{h--s)). 

By Lemma 2.1, there exists a constant, still denoted by 9, such that for 
z e Bq^, 



(3.22) w - lQ<C{ln{h-^)) 



1 X N 1 
2 
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Therefore we may choose ho small enough, so that 




^er u + M 2hl 
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which imphes ming- u > Hq. Then we finished the proof of Lemma. 

dr 

Proof of Theorem 1.1. We may assume that M — maxg- ±u, oth- 
erwise we replace u hy u — c. Then either 1 + ^ or 1 — satisfies the 
assumption of Lemma 3.4, thus Lemma 3.4 implies 

OsCq-u < (1 - ^)OsCq-u, 

which implies the C"' regularity of u near point (0, 0) by the standard iteration 
arguments. By the left invariant translation group action, we know that u is 
C° in the interior. 
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